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1. Introduction 

The pseudoscalar-meson decay constants (e.g., fn, fo,, fs and /bJ play an important role in 
extracting the CKM matrix elements (e.g., the leptonic decay width of D+ — > Z+V/ is proportional 
to /oJVcvp), which are crucial for testing the flavor sector of the standard model via the unitarity 
of CKM matrix. Experimentally, precise determination of /o and /^j^ will result from the high- 
statistics program of CLEO-c, however, the determination of /g and fs^ remains beyond the reach 
of current experiments. 

Theoretically, lattice QCD provides a solid framework to compute the masses and decay con- 
stants of pseudoscalar mesons (as well as other physical observables) nonperturbatively from the 
first principles of QCD. Thus reliable lattice QCD determinations of /b and fs^ are of fundamen- 
tal importance, in view of their experimental determinations are still lacking. Obviously, the first 
step for lattice QCD is to check whether the lattice determinations of fj) and /o^ will agree with 
those coming from the high-statistics charm program of CLEO-c. This motivated our study in 
Ref. It turns out that our predictive value of fo = 235(8)(14) MeV (posted on June 26) is in 
good agreement with the experimental result fj)+ = (223 ± 16lg) MeV announced by CLEO-c at 
Lepton-Photon Symposium (LP2005) on July 1 [|2|]. (Note that the latest value of from CLEO-c 
has been updated to /d+ = (222.6 ± 16.7+^^) MeV [|].) 

In this talk, we review our results reported in Ref. [|l|], and also present our results for the 
ratios fx/ fn and foj fu, which may be of interest from several viewpoints. 

Here we briefly outline our computations, and refer to Ref. [p (and references therein) for 
further details. First, we compute quenched quark propagators for 30 quark masses in the range 
0.03 < rtiqa < 0.80, in the framework of optimal domain- wall fermion proposed by Chiu [Q]. Then 
we determine the inverse lattice spacing a^' = 2.237(75) GeV from the pion time-con^elation 
function, with the experimental input of pion decay constant = 131 MeV. The strange quark bare 
mass nisO = 0.08 and the charm quark bare mass m^a = 0.80 ai^e fixed such that the corresponding 
masses extracted from the vector meson correlation function agree with (^(1020) and 7/ 1/^(3097) 
respectively. Then the masses and decay constants of any hadrons containing c,s, and ii{d) quarks^ 
are predictions of QCD from the first principles, with the understanding that chiral extrapolation to 
physical m,,^^ ~ nig / 25 (or equivalently m^i = 135 MeV) is required for any observables containing 
u{d) quarks. 

The observable we measure is the time-coiTclation function for pseudoscalar meson {qQ) 

Cp{t) = l^t,{y,{D, + mQ)-ly,{D, + mq)^l]^ (1.1) 

where the subscript U denotes averaging over gauge configurations. Here Cp{t) is measured for 
the following three categories: (i) Symmetric masses mg = niq, (ii) Asymmetric masses with fixed 
niQ = rus = 0.08a^^ (iii) Asymmetric masses with fixed niQ = nic = 0.80a^^ where niq is varied for 
30 masses in the range 0.03 < niqa < 0.80. Note that for mesons composed of strange and/or charm 
quarks, their masses and decay constants can be measured directly without chiral extrapolation. 

'in this paper, we work in tlie isospin limit m„ = mj. 
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The decay constant fp for a charged pseudoscalar meson P is defined by 

(0|A^(0)|P(^)>=/P^^ 

where = 75 2 is the axial- vector part of the charged weak current after a CKM matrix element 
Vqqi has been removed. Using the formula <9^A^ = {mq + mQ)qysQ, one obtains 



fp = {mq + nig) 



l(o|^75eim)l 



(1.2) 



where the pseudoscalar mass nipa and the decay amphtude z = \ (0|^75g|P(0)) | can be obtained by 
fitting the pseudoscalar time-correlation function Cp{t) to the usual formula 



2mpa 



!-mpat _^^-mpa{T-t)i 



(1.3) 



The outline of this paper is as follows. We present our results of nix and in section 2, 
niD, nto,, fo, and fo, in section 3, and the ratios fx/ fn and foJfD in section 4. In section 5, we 
summarize and conclude with some remarks. 

2. fK and mk 
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Figure 1: (a) The kaon mass niK versus the pion mass m-,i for 15 bare quark masses in the range 0.03 < 
niqa < 0. 10. The solid line is the quadratic fit. (b) The kaon decay constant fKa versus the bare quark mass 



iriqa. The solid line is the linear fit. 



We measure the time-con^elation function of kaon Cjf (f) (|0|) with rtiQ fixed at = O-OSa^^ 
while rtiq is varied for 30 masses in the range 0.03 < lUqa < 0.80. Then the data of Cjf (?) is fitted 
by the formula ( pT3| ) to extract the kaon mass niKO and the kaon decay constant /ko. 

In Fig. [T]a, the kaon mass mk is plotted versus nijc, for 15 quark masses in the range 0.03 < 
mqa<0.lO. Thedataof m/fa can be fitted by m/fa = 0.197(1) +0.255(4)(n%a) +0.389(8)(m;ta)2. 
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At the physical Umit rtiji = 135 MeV, it gives niK = 478(16) MeV, in good agreement with the 
experimental value of kaon mass (495 MeV). 

In Fig. |l|b, fxa is plotted versus bare quark mass m^a. The data is well fitted by the straight 
line //ffl = 0.068(0) + 0.1 16(1) x (m^a). At niqa = 0, it gives fK = 152(6) MeV, in agreement with 
the value fK+ = 159.8 ± 1 .4 ± 0.44 MeV complied by PDG [§1. 



3. fo, Id,, mo, and mo. 




Figure 2: (a) The D-meson mass mo versus the pion mass nij^ for 15 bare quark masses in the range 
0.03 < niqa < 0.10. The solid line is the quadratic fit. (b) The D-meson decay constant fya versus the bare 
quark mass m^a. The solid line is the linear fit. 



Now we turn to charmed pseudoscalar mesons. We measure the time-correlation function 
Cd(0 ( |1.1D with niQ fixed at = 0.80a^\ while niq is varied for 30 different masses in the range 
0.03 < niqa < 0.80. Then the data of Cz)(?) is fitted by the formula (1.3) to extract the mass mj)a 
and decay constant foa. 

In Fig. moa is plotted versus rrijca, for 15 quark masses in the range 0.03 < niqa < 0.10. 
The data of nioa can be fitted by moa = 0.816(0) +0. 101 {3){mj,a) + 0.29S{6){m^af. At m,, = 135 
MeV, it gives mjj = 1842(15) MeV, in good agreement with the mass of D meson (1865 MeV). In 
Fig. the decay constant fjja is plotted versus bare quark mass niqa. The data is well fitted by 
the straight hne /ua = 0.105(1) +0.172(1) x (m^a). Atmqa = 0, it gives /d = 235(8) MeV, which 
serves as a prediction of lattice QCD with exact chiral symmetry. 

The pseudoscalar meson of sc or sc corresponds to mga = mc-a = 0.80 and m^a = m^a = 0.08. 
Its mass and decay constant are extracted directly from the time-correlation function, which are 
plotted as the eleventh data point (counting from the smallest one) in Fig. ^. The results are 
mo,a = 0.878(2) and fo^a = 0.119(2). The mass gives mo, = 1964(5) MeV, in good agreement 
with the mass of 0^(1968). The decay constant gives fo, = 266(10) MeV, which agrees with the 
value fj)+ = 267 ± 33 MeV comphed by PDG [g]. 
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4. Ratios fK/h and foj fo 



(a) 




(b) 
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Figure 3: (a) The ratio /a-Z/tt versus r = mu/nis for 15 bare quark masses in the range 0.03 < m^a < 0.10. 
The solid line is the fit (4.1). (b) The ratio foJ fo versus mulm^. The solid line is the fit (4.2). 



At this point, it is instructive to obtain the ratios fx/ fn and foJ fo, and to see how well they 
agree the experimental values. In Fig. |3|a, the ratio ff[+ / fji is plotted versus r = niu/nis for 15 bare 
quark masses in the range 0.03 < iriqa < 0.10. The data of / fjc can be fitted by 

At physical limit m„/m.v = 1/26, it gives A = 1.17(6), in good agreement with the experiem- 
ntal value ~ 1.22 (PDG) [|]. With A = 131 MeV as the input, it gives /k^ = 153(8) MeV, in good 
agreement with the experimental value = 159.8 it 1 .4 it 0.44 MeV complied by PDG [^]. 
Similarly, in Fig. the data of fjj+ /fj)+ can be fitted by 

7.657(26) + r 

At physical limit r = 1 /26, it gives /^+ ///)+ = 1.13(2), which serves as a theoretical prediction. 

5. Summary and Concluding Remarks 

In Ref. [|], we have determined the masses and decay constants of pseudoscalar mesons K, D 
and Ds, in quenched lattice QCD with exact chiral symmetry. Our results are: 

m/f = 478 ± 16 ± 20 MeV, mo = 1 842 ± 15 ± 2 1 MeV, mo, = 1964 ± 5 ± 10 MeV, 
//f = 152±6±10MeV, /d = 235 ± 8 ± 14 MeV, /o, = 266± 10± 18 MeV, 

where in each case, the first eiTor is statistical, while the second is our crude estimate of combined 
systematic uncertainty. Note that at the time when our results were posted at arXiv:hep-ph/0506266 
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on June 26, the preliminaiy result of /o+ from CLEO-c was //)+ = 202(41) (17) MeV [Q]. So our 
value of fn was not in good agreement with the preliminary result of CLEO. However, with higher 
statistics at CLEO-c, the experimental value of fQ+ turns out to be = (222.6 ± MeV 
[^], in good agreement with our prediction of fo. Now it remains to see whether the value of fo^ 
coming from the high-statistics charm program of CLEO-c would agree with our value determined 
by lattice QCD with exact chiral symmetry. 

A salient feature of our calculation is that all quarks (no matter heavy or light) are treated on 
an equal footing, and they ai^e fully relativistic with exact chiral symmetry on the lattice. Also, 
we have not used any heavy quark approximations for the charm quark. Even though our results 
are obtained in the quenched approximation, we suspect that for lattice QCD with exact chiral 
symmetry, the quenching error in fi,, and fi,^ is less than 5%, in view of the good agreement 
between our result of fx and the experimental value. 

In closing, we note that there are several quenched/unquenched lattice QCD results for /o and 
/Dj in the literature ^, ^, as well as reported in this conference [ [Io| , For the unquenched 
staggered quark with «j = 2 + 1 [^, [I^, its prediction of fj)+ = (201 ± 3 ± 17) MeV only agrees 
with the latest CLEO experimental result ^ at 45% confidence level, and is 10% smaller than 
the experimental result. Note that this unquenched lattice QCD calculation relies on the so-called 
fourth-root trick to reduce the four tastes of each staggered quark to one taste, in computing the 
fermion determinant. Since this scheme for staggered quark has not achieved a high-precision 
prediction for fj)+ , this may be an indication of the failure of the fourth-root trick in practice. 
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